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I.  INTRODUCTION 


This  report  describes  a computer  program  that  calculates  flight 
path  data  for  an  aircraft  moving  over  the  earth.  The  program  is  called 
PROFGEN  and  was  written  in  FORTRAN.  Its  primary  intended  use  is  to 
support  simulations  that  require  a six  degree-of- freedom  trajectory 
driver. 

This  version  of  PROFGEN  evolved  from  one  written  in  1973  that 
became  obsolete  because  it  lacked  a wander-azimuth  capability  and 
employed  an  unrealistic  roll  control  mechanization.  These  shortcomings 
are  corrected  in  the  revised  version  and  several  new  features  are  added 
including  output  at  user-determined  times,  the  computation  of  attitude 
rates,  an  improved  gravity  model  and  the  ability  to  turn  through  a 
precise  angle  without  overshoot.  In  addition  the  revised  version  is 
coded  in  a modular  fashion  for  ease  of  understanding  and  change. 

This  report  will  document  PROFGEN  in  full.  Section  II  is  a 
general  description  of  PROFGEN1 s capabilities  and  limitations  that 
should  allow  the  reader  to  determine  the  program’s  applicability  to 
his  problem.  Section  III  is  a user’s  guide  that  tells  how  to  construct 
a flight  profile  with  the  available  input  parameters.  Section  IV 
develops  the  equations  that  PROFGEN  solves.  Section  V describes  the 
program  itself.  Appendix  A presents  an  example  problem  and  Appendix 
B gives  a listing  of  the  program  source  deck. 


II.  GENERAL  CHARACTERIZATION 


PROFGEN  computes  position.,  velocity,  acceleration,  attitude 
and  attitude  rate  for  an  aircraft  moving  over  the  earth.  Position 
is  given  as  (geographic)  1:  itude,  longitude  and  altitude  (see 
Figure  l).  Velocity  with  respect  to  earth  is  componentized  and 
presented  in  a local-vertical  frame  (x-y-z  in  Figure  l)  that  will 
be  called  the  navigation  frame.  Acceleration  consists  of  velocity 
rates-of-change  summed  with  Coriolis  effects  and  gravity.  Attitude 
consists  of  roll,  pitch  and  yaw,  the  Euler  angles  between  the  path 
frame  and  the  navigation  frame.  These  quantities  will  be  defined 
precisely  in  Section  IV. 

Although  the  descriptions  herein  always  refer  to  "aircraft" 
flight  paths,  PROFGEN  has  applicability  to  path  generation  for  land 
and  sea  craft  as  well.  In  general  PROFGEN  is  suited  for  simulation 
of  any  craft  under  continuous  control.  It  is  not  well  suited  to 
describing  bodies  in  free  fall  or  earth  orbit  where  mass  attraction 
is  the  primary  forcing  function. 

PROFGEN  models  a point  mass  responding  to  maneuver  commands 
specified  by  the  user.  These  maneuvers  are  available: 

• vertical  turnB  (pitch  up  or  down) 

• horizontal  turns  (yaw  left  or  right) 

• sinusoidal  heading  changes  (oscillates  left  and  right) 

• straight  flights  (great  circle  or  rhumb  line  path) 
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N-W-U  ~ Geographic  Coordinates 

x-y-z  ~ Navigation  Coordinates 

A ~ Longitude 

<f>  ~ Latitude  (geographic) 

h - Altitude 


with  exaggerated  flattening 


Figure  1 


Coordinate  Frame  Geometry 


All  horizontal-plane  maneuvers  are  executed  in  a coordinated  fashion. 

This  simply  means  that  the  aircraft  is  rolled  to  an  angle  where  the 

vector  sum  of  the  centrifugal  turning  force  and  the  force  of  "gravity" 

2 

(32.2  ft /sec  ) acts  perpendicular  to  the  wings.  Only  one  type  of 
maneuver  may  he  executing  at  any  given  time  but  it  can  commence  from 
any  aircraft  attitude.  For  example,  the  aircraft  may  go  into  a left 
turn  while  in  a dive. 

In  addition  to  the  four  basic  maneuvers,  the  user  also  has 
control  of  path  acceleration  by  which  the  aircraft  can  be  forced  to 
change  speeds.  Path  acceleration  may  be  superimposed  over  any 
maneuver.  This  would  allow,  for  example,  an  accelerated  diving  turn. 

PROFGEN  is  used  to  create  an  extended  flight  profile  by  stringing 
together  a sequence  of  maneuvers  chosen  from  the  basic  four.  The  user 
specifies  how  long  each  maneuver  shall  last  and  thereby  divides  the 
profile  into  flight  segments . Up  to  fifty  flight  segments  may  be 
strung  together  to  produce  a varied  total  profile.  The  final  values 
of  the  variables  in  each  segment  are  passed  along  as  the  initial  values 
for  the  start  of  the  next  segment  thereby  creating  uninterrupted  time 
histories  for  all  output  variables. 

The  program  allows  step  changes  to  occur  in  displacement  accel- 
eration and  in  rotational  velocity.  This  produces  continuous  time 
histories  for  displacement  velocity  and  rotational  position  (roll, 
pitch,  yaw)  but  results  in  infinite  Jerk  (rate-of-change  of  displace- 
ment acceleration)  and  infinite  rotational  acceleration. 
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Acceleration,  velocity  and  position  are  related  instantaneously 
by  integration  and  differentiation  to  within  the  accuracy  of  the 
Kutta-Merson  numerical  integrator.  Every  effort  has  been  made  to 
configure  this  integrator  to  produce  an  accurate  result  so  that  the 
output  variables  form  a self-consistent  set.  Thus  the  integrator  is 
fifth  order  and  can  adjust  its  step  size  automatically  to  control 
the  growth  of  errors.  To  illustrate,  a great  circle  path  from  Dayton  to 
Moscow  accumulated  less  than  15  feet  of  error  over  its  5000  mile  distance. 

PROFGEN  is  limited  in  its  capability  to  simulate  intricate 
fighter  maneuvers.  This  arises  in  part  because  PROFGEN  forces  the 
aircraft  body  and  path  frames  to  be  coincident  and  thereby  loses 
the  ability  to  simulate  slipping  or  crabbing  motion.  Thus,  for 
example,  one  could  not  simulate  a fighter  aircraft  doing  a barrel 
roll  or  an  Immelmann,  On  the  other  hand,  one  ,ould  simulate  a 
complete  loop  of  arbitrary  radius  since  severity  of  maneuver  is 
not  restricted.  In  general,  PROFGEN  can  simulate  any  maneuver 
possible  with  a bomber  or  cargo  aircraft. 

The  earth  is  modeled  as  a perfect  ellipsoid  naving  values  for 
eccentricity,  semimajor  axis  length,  spin  velocity  and  gravitational 
constant  equal  to  those  of  the  DOD  World  Geodetic  System  1972 
(Ref.  l).  Earth’s  gravity  is  modeled  as  a function  of  latitude 
and  altitude,  having  both  radial  and  level  components.  This  model 
is  not  overly  precise  (probably  no  better  than  25  micro  gees) 
and  may  need  revision  for  some  applications. 
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PROFGEN  compiles  and  executes  in  less  them  60,000  words  of 
CDC  CYBER-7I*  memory.  It  uses  only  single  precision  variables  and 
slLI  source  code  is  FORTRAN.  The  program  takes  six  seconds  of 
central  processor  time  to  compile.  The  ratio  of  simulated  time 
to  execution  time  improves  as  problem  dynamics  become  less  severe, 
reaching  20267  : 1 for  straight  flight  segments  but  falling  to 
h : 1 for  a 10  gee  horizontal  turn. 
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III.  USER'S  GUIDE 


This  section  defines  the  input  data  that  the  user  supplies  to 
run  PROFGEN.  The  input  data  specifies 

• initial  conditions 

• maneuver  characteristics 

• integrator  control 

• output  control 

All  data  is  entered  under  a NAMELIST  format  that  permits  the 
entry  of  character  strings.  A character  string  is  a parameter  name 
followed  hy  its  values  written  in  the  user's  choice  of  format 
specification.  The  use  of  NAMELIST  on  the  CDC  CYBER-7^  will  be 
illustrated  In  Figures  2 and  3. 

Two  NAMELIST  input  data  lists  are  used,  PRDATA  and  PASDATA.  The 
PRDATA  (Problem  Data)  list  contains  15  parameters  that  remain  fixed 
for  the  entire  run.  These  parameters  specify  all  initial  conditions 
and  control  output. 

The  PASDATA  (Path  £egment  Data)  list  contains  13  parameters  that 
remain  fixed  only  for  the  length  of  a segment.  These  parameters  specify 
and  describe  each  maneuver,  control  the  numerical  integrator,  and 
control  the  output  frequency. 

3.1  FRDATA  Input 

Fifteen  parameters  are  entered  through  the  PRDATA  list.  Failure 
to  specify  any  one  of  these  results  in  program  termination.  All 


parameters  sire  single  precision  and  all  must  be  entered  in  units  of 
feet,  seconds  and/or  degrees.  The  following  format  will  be  used 
to  describe  input  peu*ameters  throughout  this  section  euid  the  next. 


Parameter 


(Type) 


Units  (If  Any) 


IPROB  (Integer) 

The  problem  identification  number.  It  is  set  by  the  user 
for  identification  purposes  only. 


NSEGT  (Integer) 

The  total  number  of  path  segments  required  to  complete  the 
entire  problem.  This  number  may  not  exceed  50  as  the  program 
is  now  configured. 


LLMECH  (Integer) 

The  local-level  azimuth  .single  mechanization  index. 
See  Section  and  Table  2. 

LLMECH  Azimuth  Mechanization 

1 Alpha  Wander 

2 Const  suit  Alpha 

3 Unipolar 

U Free  Azimuth 


TSTART  (Real)  seconds 

The  initial  time.  It  is  used  to  begin  the  problem  at  any 
desired  point.  It  may  be  negative. 
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VTO  (Real)  feet  per  second 

The  initial  magnitude  of  total  velocity  with-respect-to  the 
earth.  VTO  must  be  non-negative. 


PHEADO  (Real)  degrees 

The  initial  heading  angle  of  the  path  coordinate  frame. 
It  is  specified  as  positive  clockwise  from  North.  Its 
range  is  the  closed  interval  [-180.,  +180.]. 


PPITCHO  (Real)  degrees 

The  initial  pitch  angle  of  the  path  coordinate  frame.  It 
is  specified  as  positive  in  the  upward  direction.  The  path 
frame  is  level  when  the  pitch  angle  is  zero.  Its  range  iB 
[-90.,  +90.]. 


ALFAO  (Real)  degrees 

The  initial  alpha  angle.  Alpha  is  the  navigation  frame 
heading  angle  and  is  specified  positive  counterclockwise 
from  North.  Its  range  is  [-180.,  +l80.j. 


LATO  (Real)  degrees 

The  initial  geographic  latitude.  Its  range  is  the  open 
interval  (-90.,  +90.).  Since  the  program  falters  when 
trying  to  compute  at  exactly  90  degrees,  these  two  extreme 
points  must  be  avoided. 
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LONO 


(Real) 


degrees 


AXiTO 


IPRNT 


The  initial  longitude.  It  has  no  effect  on  the  problem’ a 
dynamics  but  is  necessary  to  establish  a reference  point 
for  the  calculation  of  current  position.  Its  range  is 
[-180.,  +180.3. 


(Real)  feet 

The  initial  altitude  above  the  reference  ellipsoid. 
ALTO  may  be  negative. 


(Integer) 

Print  control  index  having  control,  in  part,  over  what  is 
written  on  TAPE6.  This  tape  is  considered  to  be  printed 
output.  All  TAPE6  output  is  formatted. 

IPRNT  Action 

1 Output  on  TAPE6  at  time-intervals  specified 

by  DTO  (a  PAS DATA  parameter) 

1*1  Output  at  DTO  intervals  is  turned  off. 

Regardless  of  the  state  of  IPRNT,  the  following  output  also 
appears  on  TAFE6: 

. date  and  time 

• input  data  from  PRDATA  and  PASDATA  lists 

• variable  values  at  start  of  each  segment  and  at  t-final 

• error  messages 

. post-run  assessment  of  numerical  integrator  performance 
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IRITE  (Integer) 

Write  control  index.  This  output  is  written  on  TAPE3  and 
is  designed  for  compact  storage  of  data' for  subsequent  use 
by  another  program.  All  TAPE3  output  is  unformatted. 


IRITE  Action 

1 Output  on  TAPE3  consisting  of  date,  time,  input 

data  and  variable  values  beginning  at  TSTART 
and  continuing  at  DTO  intervals. 

*1  No  output  on  TAPE3. 


IPLOT  (Integer) 

Plot  control  index.  This  output  is  on  PLFILE  for  post-run 
graphing  using  DISSPLA,  a CALCOMP  plot  library. 

IPLOT  Action 

1 Program  plots  five  graphs,  latitude  vs.  longitude 

and  time  histories  of  altitude,  roll,  pitch  and 

yaw.  Up  to  501  points  are  plotted  in  each  graph, 
£he  first  being  at  TSTART  and  all  thereafter 
at  ITO  intervals. 

'1  No  plotted  output. 


ROLRATE  (Real)  degrees  per  second 

Nominal  aircraf*  roll  rate.  When  the  aircraft  must  bank  to 
execute  b coordinated  horizontal  turn,  it  rolls  to  the  proper 
bank  angJe  at  a rate  of  ROLRATF.  In  sine-heading-change 
maneuvers,  ROLRATE  serves  as  the  limiting  value  for  the 
derivative  of  roll.  ROLRATE  must  be  positive. 
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Figure  2 is  a sample  of  a PRDATA  card  input  set.  Note  that  the 
data  items  may  he  listed  in  ary  order  so  long  as  they  all  appear 
between  the  beginning  identifier*  $ PRDATA,  and  the  ending  identifier, 
$. 


3.2  PASDATA  Input 

Thirteen  parameters  having  up  to  50  values  each  are  entered 
through  the  PASDATA  list.  Each  parameter  is  dimensioned  in  the 
program  as  a 50  element  array,  the  number  50  corresponding  to  the 
maximum  number  of  segments  allowed.  Each  parameter  value  must  be 
assigned  to  the  array  element  corresponding  to  its  segment  number; 
for  example,  if  the  output  spacing  in  the  sixth  segment  is  to  be 
25  seconds,  one  would  input  DT0(6)  * 25.  Fach  parameter  name  in 
the  list  that  follows  has  the  argument  i appended  to  it  to  indicate 
its  dependence  on  segment  i,  1 — i — 50. 

Six  of  the  PASDATA  parameters  (TURN,  NPATH,  PACC,  TACC,  HEAD, 
PITCH)  describe  the  maneuver  and  four  (MODE,  ERROR,  HMAX,  HMIN)  are 
associated  with  numerical  integration.  The  other  three  control  output 
frequency  (PTO),  set  segment  length  (SEGLNT) , and  control  initial 
conditions  (RESTART),  Each  parameter  has  a default  option  that  is 
invoked  in  lieu  of  input  data.  The  default  saves  the  user  the 
trouble  of  specifying  values  that  often  recur.  All  parameters  are 
single  precision  and  all  must  be  entered  in  units  of  feet,  seconds, 
gees  (l  gee  = 32.2  ft/sec.  ) and/or  degrees. 


JPROATA  IPR0B*65Q, 


NSEGT*i?» 

LLMECH*2» 

TST  ART=0.  » 
VTO»1000 • • 

PHE  AOOS180* * 
PPI TCHO*Q • « 
ALFAO=<*5  • t 
ALT 0=30000.* 
LArO*39.t 
L0N0=-84. » 
ROLRATE«250.  * 
IPRNT=l, 
IRITE*0* 

IPLOT *1S 


Figure  2 - Sample  of  PREATA  Input 


Parameter 


vTypj)  Units  (if  Any) 


SEOLNT(i)  (Real)  seconds 

XU 

The  time  interval  of  the  i segment.  SEGLNT(i)  cam  be 
any  non-negative  number,  including  zero.  The  program 
remains  in  segment  i until  exactly  SEGLNT(i)  seconds  have 
been  simulated.  The  default  value  is  zero  seconds. 


RESTART(i) 


( Integer ) 


The  index  number  for  control  of  the  initial  conditions  at 
the  beginning  of  each  segment. 

RESTRART(i)  Action 

1 All  variables  in  the  state  vector  are  reset 

to  the  conditions  that  existed  at  TSTART, 
namely  those  in  PRDATA.  RESTART  * 1 is 
useful  when  one  wishes  to  produce  a reference 
flight,  and  a variation  of  that  flight,  all 
in  one  run. 

J*1  The  variable  values  at  the  beginning  of 

segment  i equal  those  at  the  end  of 
segment  i-1. 

The  default  value  is  zero,  no  reset  performed. 


TURN(i) 

The  index  number  for  the  type  of  maneuver  to  be  used. 
TURN(i)  Action 

1 vertical  turn 

2 horizontal  turn 

3 sinusoidal  heading  change 

4 


straight  flight 


All  maneuvers  begin  at  the  start  of  a segment.  Vertical  and 
horizontal  turns  are  complete  when  a specified  turn  angle  is 
reached.  If  specified  angle  is  reached  and  time  remains  in 
the  segment,  PROFOEN  reverts  to  a straight  flight  mode 
(TURN  * It)  for  the  remaining  seconds  of  the  segment.  If 
TURN(i)  is  3,  a '’sinusoidal”  path  (oscillatory  yawing  motion 
in  the  horizontal  plane)  is  flown  for  SEGLNT(i)  seconds. 

For  sine  maneuvers,  the  user  must  select  a segment  length 
that  is  a multiple  of  TpA  where  Tp  is  the  period  of  the 
sinusoid.  If  TURN(i)  is  4,  a straight-flight  segment  will 
be  flown  over  a nominal  path  determined  by  the  value  of 
NPATH(i)  for  SEGLNT(i)  seconds.  Section  3A  discusses 
these  maneuver  characteristics  more  fully.  The  default 
value  is  U,  straight  flight. 


NPATH(i)  (Integer) 

The  index  number  for  the  nominal  path. 

NFATH(i)  Action 

1 Great  cricle  path 

2 Rhumb  line  path 

When  a rhumb  line  path  is  chosen,  the  aircraft  maintains  a 
constant  heading  angle  during  straight  flight  periods.  When 
a great  circle  path  is  chosen, the  aircraft  flies  in  a fixed 
plane  during  straight  flight  periods.  The  aircraft  maintains 
this  fixed-plane  flight  over  the  ellipsoidal  earth,  even  when 
altitude  changes,  by  correcting  heading  continuously.  When 
not  in  straight  flight  (i.e.  TURN  = 1,  2 or  3),  the  rhumb 
line  or  great  circle  actions  are  superimposed  on  the  chosen 
maneuver.  The  default  value  is  2,  rhumb  line  path. 


i 

l 


•i 


. 

3 


PACC(i)  (Real)  gees 

The  signed  value  of  the  constant  acceleration  along  the 
velocity  vector,  i.e.  along  the  path  x-axis.  The  program  2 
converts  PACC(i)  in  gees  to  path  acceleration  in  feet/second 
by  multiplying  by  32.2.  PACC(i)  may  be  assigned  any  real 
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value;  it  remains  that  value  for  the  entire  segment 
regardless  of  maneuver  specification.  Positive  (negative) 
values  cause  the  aircraft  to  gain  (lose)  total  speed. 

Since  all  active  maneuvers  (TURN  =1,  2 or  3)  require  a 
division  by  total  speed  (VT^  to  compute  acceleration,  the 
user  must  assign  PACC(i)  so  VT  is  never  zero  during  the 
actual  turning  portion  of  such  maneuvers.  PACC(i)  may 
force  VT  to  zero  anytime  during  a straight  flight  segment. 
The  default  is  zero  gees. 


TACC(i)  (Real)  gees 

« 

The  magnitude  of  the  maximum  centrifugal  acceleration 
during  either  a vertical  or  horizontal  turn.  The  program 
converts  TACC(i)  in  gees  to  acceleration  in  feet/second^ 
by  multiplying  by  32.2.  TACC(i)  must  be  positive  for 
vertical  and  horizontal  turns t The  default  value  is  zero 
gees. 


HEAD(i)  (Real)  degrees 


! 


HEAD(i)  has  two  uses.  j 

\ 

For  horizontal  turns,  HEAD(i)  is  the  desired  change  in 

heading  angle.  Other  factors  permitting  (SEGLNT,  TACC,  | 

ROLRATE,  PACC,  VT)  this  turn  angle  will  be  executed  f 

accurately.  The  magnitude  of  HEAD(i)  may  be  greater 

than  360  degrees.  A positive  (negative)  HEAD(i)  forces  ] 

a right  (left)  turn.  1 

j 

For  sine  maneuvers  HEAD(i)  is  the  maximum  variation  of  j 

the  heading  angle  and  its  absolute  value  must  be  less  ; 

than  90  degrees.  A positive  (negative)  HEAD(i)  forces  \ 

the  sine  maneuver's  ground  track  to  lie  right  (left)  * 

of  the  initial  ground  track.  The  default  value  is  ) 

zero  degrees . ; 


PITCH(i) 


(Real) 


degreeB  or  deg/sec 


PITCH(i)  has  two  uses. 


For  vertical  turns  PlTCH(i)  is  the  desired  change  in  pitch 
angle  in  degrees.  Other  factors  permitting  (SEGENT,  TACC, 

PACC,  V ),  this  value  will  "be  achieved  precisely.  PIrj.CH(i) 
may  exceed  90  degrees.  A positive  (negative)  PITCH(i) 
forces  the  pitch  angle  to  increase  (decrease). 

For  sine  maneuvers,  PITCH(i)  is  the  frequency  of  the  sinusoidal 
rate  of  change  of  heading  in  degrees  per  second.  It  must 
be  non-zero.  The  sign  of  PITCH(i)  has  no  effect  on  the 
sine  maneuver.  The  default  value  is  zero  in  degrees  or 
degrees  per  second,  as  the  case  may  be. 


DTO(i) 


(Real) 


seconds 


The  time  interval  between  required  output  times.  DTO(i) 
is  referenced  to  zero  seconds;  e.g,,  if  DTO(i)  = 6,  output 
would  be  available  at  T » (•••,  -12,  -6,  0,  6,  12, • • • ). 
DTO(i)  must  be  positive.  DTO(i)  controls  output  frequency 
for  printing,  writing  and  plotting  (see  IPRNT,  IRITE,  IPLOT). 
Careful  sizing  of  DTO(i)  is  a necessity,  especially  when 
two  or  three  output  modes  are  used  simultaneously.  The 
default  value  is  100  million  seconds  corresponding  to  no 
output  at  all. 


MODE(i)  (Integer) 

The  index  for  control  of  step  size  In  the  numerical 
integration  routine. 


MODE(i) 

0 

1 


Action 

Fixed  step-size  integration. 
Variable  step-size  integration. 


The  step  size  is  HMIN(i)  when  fixed  step-size  integration 
is  used.  A fifth  order  numerical  integration  is  performed. 


IT 


v 


With  the  variable  step-size  mode,  the  program  begins  the 
integration  with  a step  size  of  HMIN(i).  The  numerical 
integrator  adjusts  the  step  size  upwards  from  there  while 
keeping  the  within-step  error  below  the  value  specified 
in  ERROR(i).  If  problem  dynamics  are  mild,  the  step  size 
can  grow  very  large,  limited  finally  by  HMAX(i).  If  problem 
dynamics  are  severe,  the  minimum  step  size  may  not  be 
adequately  small  to  satisfy  the  error  criterion  in  which 
case  an  error  message  is  printed. 

In  summary  both  integration  modes  perform  fifth  order 
numerical  integrations  but  MODE  = 1 adjusts  step  size 
automatically  to  conform  to  an  error  criterion.  The 
default  value  is  variable  step-size  integration. 


ERROR ( i ) (Real) 

The  allowable  within-step  integratign  error.  It  must  be 
positive.  The  default  value  is  10~  , a value  that  has 
proven  satisfactory  during  program  development. 


HMAX(i)  (Real)  seconds 

The  maximum  step  size  when  variable  step-size  integration 
is  used.  It  must  be  positive.  The  default  value  is 
10,000  seconds. 


HMIN(i)  (Real)  seconds 

The  minimum  step  size  when  variable  step-size  integration 
is  used.  With  fixed  step-size  integration,  HMIN(i)  is 
the  size  of  each  step.  It  must  be  positive.  The  default 
value  is  one  second. 


Table  1 shows  the  relationship  of  TURN,  TACC,  HEAD  and  PITCH. 
Figure  3 is  a sample  of  a PASDATA  card  input  set.  Note  that  some 
parameters  are  not  specified  because  the  desired  values  agreed  with 
the  default  option.  Also  note  the  capability  to  specify  repeated 
values  using  a repetition  factor. 
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figure  3 - Sample  of  PASDATA  Input 


TABLE  1 DEFINITION  OF  TURN  P 


♦See  text  for  units 


3*3  Program  Limitations  (What  Happens  If  . . . ) 

PROFGEN  will  not  begin  profile  generation  until  each  parameter 
lies  within  its  permitted  range  as  specified  in  3.1  and  3.2.  Subroutine 
VALDATA  range-checks  NSEGT,  LLMECH,  VTO,  PHEADO,  PPITCHO,  ALFAO,  LATO, 
LONO,  ROLRATE,  SEGLNT,  TURN,  NPATH,  TACC,  HEAD,  PITCH,  DTO,  MODE, 

ERROR,  HMAX  and  HMIN.  A message  is  printed  for  each  range-check 
that  fails  and  the  program  is  terminated. 

Error  messages  can  also  occur  during  profile  generation  (i.e. 
after  TSTART).  One  such  mid-run  message  occurs  if  and  when  the 
integrator  reduces  step  size  to  HMIN  and  is  still  not  able  to  satisfy 
the  error  criterion  (ERROR).  In  such  cases  this  message  is  printed: 

THE  INTEGRATION  ERROR  EXCEEDS  ITS  ALLOWED  VALUE 
When  this  occurs  PROFGEN  is  designed  to  continue  to  run,  doing  the 
best  it  can  with  HMIN.  The  value  of  the  result  is  questionable, 
however,  and  the  best  advice  is  to  scrap  the  output,  reduce  HMIN 
by  at  least  a factor  of  ten,  and  rerun  the  program. 

Another  mid-run  error  message  occurs  if  and  when  the  product  of 
computed  roll  rate  and  minimum  step  size  would  produce  a roll  bank 
angle  in  excess  of  90  degrees.  Since  the  aircraft  must  bank  to 
execute  either  a horizontal  turn  or  a sine  maneuver,  excessive  roll 
angles  could  occur  in  either  type  of  maneuver.  PROFGEN  avoids  this 
problem  in  a horizontal  turn  but  succumbs  to  it  in  a sine  maneuver; 
prior  to  each  sine  maneuver  the  program  checks  for  the  problem  and, 
if  it  exists,  prints  the  following  warning  message  and  then  terminates 
execution. 
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CHKSHC  MESSAGE  - THE  PRODUCT  OF  COMPUTED 


ROLL  RATE  AND  MINIMUM  STEP  SIZE  EXCEEDS 
90  DEGREES.  BANK  ANGLES  IN  EXCESS  OF 
90  DEGREES  ARE  NOT  ALLOWED.  PROGRAM 
TERMINATED. 

Again  the  solution  is  to  reduce  HMIN  for  that  segment. 

Another  mid-run  message  occurs  if  and  when  the  cosine  of  pitch 
is  exactly  zero.  This  would  happen,  of  course,  if  pitch  magnitude 
were  exactly  JI/2  radians  (90  degrees).  At  90  degrees,  the  algorithm 

for  computing  yaw  rate  and  roll  rate  would  make  both  of  these 

\ 

quantities  infinite.  PROFGEN  recognizes  the  situation  and  prints 
the  following  warning  message  from  subroutine  ETADOT. 

ROLL  AND  YAW  RATES  ARE  UNDEFINED 
WHEN  PITCH  IS  90  DEGREES.  THUS 
ALL  RATES  HAVE  BEEN  TEMPORARILY 
ZEROED. 

No  divisions  by  zero  are  attempted  so  the  program  continues  to 
execute.  In  short  PROFGEN  handles  a pitch  angle  of  1 90  degrees 
by  avoiding  the  fatal  rate  computations. 

If  latitude  becomes  1 90  degrees,  PROFGEN  attempts  a division  by 
zero  in  LAMD0T  and  suffers  a fatal  error  in  which  the  CDC  operating 
system  kicks  the  program  off  the  machine.  Similar  zero-division 
failures  occur  when  one  attempts  a horizontal  plane  maneuver 
(horizontal  turn  or  sine  maneuver)  with  horizontal  velocity  equal 
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zero,  or  when  a vertical  turn  is  attempted  with  total  velocity  equal 
zero,  or  when  the  aircraft  is  flovn  into  the  earth’s  center.  Other 
zero-division  situations  would  he  even  rarer  than  these  and  are  not 
worth  mentioning. 

3.4  What  to  Expect  from  Each  Maneuver 

This  section  describes  each  maneuver  in  depth  to  see  what  it 
does  and  how  it  does  it.  These  descriptions  form  the  basis  for 
the  development  of  the  control  equations  in  Section  4.3. 

3.^.1  Verticail  Turn 

A vertical  turn  is  a pitch-up  or  pitch-down  maneuver  that  takes 

place  in  a vertical  plane.  As  with  all  maneuvers,  vertical  turns 

begin  executing  at  the  start  of  a segment  (TI).  Pitch  angle  advances, 

at  a rate  controlled  by  TACC  and  aircraft  speed,  until  the  time  in 

the  segment  runs  out  at  TF  or  until  the  change-in-pitch  reaches  PITCH 

degrees  at  TDONE,  whichever  time  comes  first.  Altitude,  pitch  and 

acceleration  curves  for  two  vertical  turns  are  shown  in  Figure  4. 

Let  a^  represent  turn  acceleration  normal  to  the  flight  path. 

PROFGEN  holds  a ( “TACC  fps2  ) constant  while  pitch  advances.  Since 
n 

V* 

a,  « — - Ve  o> 

the  turn's  radius  of  curvature,  r,  and  its  advancement  rate,  0, 
are  also  constant  as  long  as  total  speed,  V,  remains  fixed. 
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Altitude 


Turning  action  is  enabled  by  switching  6 on  at  TI  and  then  off 
at  min  (TF,  TDONE).  This  produces  a pitch-rate  discontinuity  at 
min  (TF,  TDONE)  that  the  numerical  integrator,  KUTMER,  cannot 
handle.  PROFGEN  solves  the  problem  by  splitting  the  segment  into 
two  piecea  one  from  TI  to  TDONE  and  the  other  from  TDONE  to  TF. 

(If  TDONE  > TF,  only  one  piece  is  necessary,  viz.  TI  to  TF.) 

KUTMER  integrates  the  two  disjoint  pieces  separately  and  thereby 
avoids  a time  step  that  would  span  the  pitch-rate  discontinuity. 

The  switching  action  on  8 may  be  observed  in  the  program's 
pitch-rate  output  which  is  a non-zero  constant  while  pitch  is 
advancing  and  zero  thereafter.  Vertical  plane  maneuvers  induce 
no  rolling  or  yawing  motion. 

TDONE  is  computed  in  subroutine  TSETUP1  before  segment 
integration  begins.  The  computation  for  TDONE  assumes  two  things: 

• turn  acceleration  is  constant 
. total  speed  does  not  drop  to  zero 
The  first  assumption  is  guaranteed  by  the  program's  construction. 

The  user  must  guarantee  the  second  assumption  by  choosing  PACC  so 
total  speed  will  remain  positive.  When  these  assumptions  hold 
the  aircraft's  PITCH  angle  will  advance  exactly  PITCH  degrees 
in  the  interval  TI  to  TDONE  as  illustrated  in  Example  1 of  Figure  U, 
If  TDONE  exceeds  TI,  the  change-in-pitch  will  fall  short  of  PITCH  as 
illustrated  in  Example  2 of  Figure  1*. 

The  minimum  time  required  to  complete  a vertical  turn  through 
an  arbitrary  pitch  angle  A0  is  as  follows: 


At 


, N4-o 

a* 


(2) 


where  At  * time  required  to  pitch  through  A6  radians  ( >0) 

Vq  ■ total  speed  at  TI  (>o) 

A0  ■ turn  angle  ■ |PITCH|  ( >o ) 

an  * normal  turning  acceleration  ■ TACC  (>o) 

• 

Vq  ■ tangential  acceleration  ■ PACC 

A derivation  of  this  result  is  given  in  Section  4.3.2.  Equation  (2) 
is  useful  for  computing  flight  time  in  a pitch  maneuver. 


3.4.2  Horizontal  Turn 

In  a horizontal  turn  the  aircraft  heading  swings  left  or  right 
to  force  the  aircraft  to  follow  a pseudo-circular  path  over  the 
ground.  Such  a turn  can  he  performed  in  any  pitch  attitude  except 
± 90  degrees.  Horizontal  turns  are  always  performed  in  coordinated 
fashion.  (Coordinated  turns  are  also  termed  symmetric.)  A coordinated 
turn  is  one  in  which  the  aircraft  roll  (bank)  angle  ic  controlled 
so  that  the  vector  sum  of  the  horizontal  turning  force  and  the 
vertical  force  of  "gravity"  (defined  for  this  purpose  as  32.2  ft/sec  ) 
acts  perpendicular  to  the  wingB,  For  example,  in  a level  one-gee 
turn  to  the  pilot's  right,  the  aircraft  rolls  about  its  long  axis 
to  a bank  angle  of  45  degrees,  right  wing  down.  Because  heading  and 
roll  must  both  be  controlled,  the  software  implementation  for  the 
horizontal  turn  is  more  complex  than  that  for  the  vertical  turn. 
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As  was  true  with  pitch  in  the  vertical  turn,  heading  advances 


'A 

I 
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in  the  horizontal  turn  until  the  tine  in  the  segment  runs  out  at 
TF  or  until  the  change-in-heading  reaches  HEAD  degrees  at  TDONE, 
whichever  tine  comes  first.  Another  way  to  say  this  is  that  the 
aircraft  turns  in  the  time  interval  between  TI  and  min  (TF,  TDONE). 
During  this  turning  interval,  while  heading  advances  continuously, 
roll  also  goes  through  its  own  set  of  gyrations  in  order  to 
implement  a coordinated  turn.  Representative  roll  curves  are  shown 
in  Figure  5- 

Note  that  roll  always  begins  and  ends  at  zero  and  remains  in 

the  interval  (-90°,  +90°).  Also  note  that  when  roll  changes,  it 

does  so  at  the  constant  rate,  ROLRATE. 

In  contrast  to  the  vertical  turn  where  a was  constant,  a for 

n n 

the  horizontal  turn  follows  a curve  similar  in  shape  to  the  roll 
curves  from  Figure  5.  a^  is  given  by 

an  ct)  - 3z.z  x**  At  ( nk  ct))  (3> 

where  n is  (constant)  pitch  and  ri  is  roll.  Note  that,  since  n 
y x x 

varies  with  time,  a does  also  thereby  producing  a path  with  a 

n 

variable  radius  of  curvature.  (The  radius  of  curvature  is  infinite 
at  the  two  ends  of  the  turn  and  reaches  a minimum  when  bank  angle 
peaks.)  Lat-long,  yaw,  roll  and  acceleration  curves  for  two 
horizontal  turns  are  shown  in  Figure  6. 
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It  is  apparent  from  Figure  5 that  roll  rate  has  from  one  to 
three  points  of  discontinuity  within  the  segment  - one  in  Case  D, 
two  in  B and  C and  three  in  A.  Again,  the  numerical  integration 
problem  that  this  presents  is  handled  by  piecewise  integration  as 
explained  in  Section  3.^.1. 

Before  integration  begins,  time  points  TOFF,  TON  and  TDONE 
(defined  in  Fig.  5)  are  computed  in  subroutine  TSETUP2.  The 
condition  on  TDONE  is  that  heading  at  TDONE  should  be  different 
from  heading  at  TI  by  HEAD  degrees.  To  compute  TDONE,  TSETUP2  must 
account  for  variations  in  both  acceleration  (an(t))  and  speed.  The 
exf  t equations  for  doing  this  are  very  non-linear  and  have  been 
approximated  in  PROFGEN  as  quadratics  in  TDONE.  If  TSETUP2  finds 
TDONE  is  larger  than  TF  it  makes  TDONE  equal  to  TF  to  keep  the 
turn  within  the  time  limit  of  the  segment.  Once  TDONE  is  known, 

TOFF  and  TON  are  easily  computed  based  on  max  roll  angle  and 
ROLRATE.  As  in  the  vertical  turn,  PROFGEN  assumes  that  speed 
remains  positive  throughout  the  turn  segment , a condition  that 
the  user  must,  guarantee. 

The  following  equation  is  an  approximate  expression  for  the  time 
required  to  complete  a turn  through  Alp  radians. 


where 


At  * time  required  to  turn  Alp  radians  (>o) 

Vq  = total  speed  at  TI  (>o) 

Aip  - turn  angle  = |HEAD|  (>o) 

a^  * normal  turning  acceleration  = TACC  (>o) 

• 

Vq  = tangential  acceleration  = PACC 

2(T0FF-TI)  = time  required  to  roll  into  and  out  of  turn 


= 2 tan-1  t ^ /ROLRATE 

132.2  cos(TX^)  )/ 

Thi;i  equation  is  approximately  correct  for  a turn  that  rolls  quickly 
to  its  maximum  hank  angle,  holds  that  angle  for  awhile  and  then  rolls 
quickly  hack  to  zero  (Case  A in  Figure  5 ).  The  error  in  this  equation 
grows  large  as  A\p  and  ROLRATE  grow  smaller  and  as  PACC  and  TACC  grow 
larger. 

3.^.3  Sine  Maneuver 

In  a sine  maneuver  the  aircraft  follows  a ground  path  like  that  of 
Figure  7a.  This  path  results  when  ground  heading,  ^(t),  is  controlled 
hy  the  equation 


V(t)  - f +A 


o it  < Tp/z 


- A > Tpfe  4 t < 7p  s 


where  A is  maximum  heading  variation  (HEAD)  and  w is  oscillation 
frequency  (PITCH). 
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Repeated  cyles  of  7a  are  shown  in  7b  and  are  produced  by  simply 
iterating  the  above  equation  to  yield  a longer  maneuver  similar  to 
jinking.  Note  that  neither  7a  or  7b  are  properly  scaled. 

A sine  maneuver  may  execute  in  any  pitch  attitude  except  ±90 
degrees  and  is  always  performed  in  coordinated  fashion.  Again,heading 
and  roll  must  both  be  controlled  but  the  governing  equation  is  the 
one  for  heading  given  above.  The  companion  equation  for  roll  that 
produces  coordinated  maneuvers  is 

y ifinr  ***  tv 

where  V is  total  speed.  Since  r\  has  no  discontinuities, the  numerical 

i * 

integration  can  proceed  uninterrupted  the  sine  maneuver  thereby 

avoids  complex  event-time  calculations  like  those  for  a horizontal  turn. 

Figure  8 shows  ground  track,  roll  and  heading  curves  (to  scale)  for 

sine  maneuver  where  A is  -20°,  T is  10  seconds,  V is  1000  fps  and 

P 

SEGLNT  is  12.5  seconds.  Note  that  roll  passes  through  zero  at 
multiples  of  T A seconds  so  that  the  aircrafts  wings  are  level  when 
the  segment  is  finished  at  12.5  seconds. 

3A.1*  Straight  Flight 

Complete  straight-flight  segments  occur  when  TURN  is  It  and  partial 
segments  occur  anytime  a vertical  or  horizontal  turn  has  reached  its 
max  turn  angle  with  time  remaining  in  the  segment.  Neither  roll  nor 


gure  r 
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Example  of  Sine  Maneuver 
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pitch  vary  in  straight  flight  segments  and  heading  is  governed  by  the 
users  choice  of  nominal  path  (NPATH).  Heading  is  constant  over  a 
rhumb  line  path  whereas,  for  a great  circle  path,  heading  must  vary 
to  keep  the  aircraft  in  the  great  circle  plane.  Rhumb  line  flights 
that  continue  long  enough  spiral  in  on  one  of  the  earth’s  poles  and 
end  up  causing  a division-by-zero  failure. 

Total  speed,  which  had  to  remain  positive  during  turning  maneuvers 
may  be  zero  in  straight  flight  segments.  At  such  times,  aircraft 
position  is  fixed  and  attitude  is  that  which  existed  just  prior  to 
speed  becoming  zero. 

To  aid  the  user  in  constructing  straight  flight  segments  between 
locations  over  the  earth,  a program  called  HEADING  has  been  written. 

In  response  to  user  inputs  of  lat,  Ion  and  altitude  at  origin  and 
destination,  HEADING  computes  the  heading  angle  at  origin  neei  i to 
reach  destination  over  a great  circle  path.  HEADING  also  computes  the 
great  circle  distance  from  origin  to  destination.  HEADING  is  a 
double  precision  FORTRAN  program  that  can  be  made  available  to 
interested  users. 
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IV.  ANALYTICAL  DEVELOPMENT 


This  section  develops  the  equations  that  govern  the  trajectory 
of  an  aircraft  under  continuous  control  in  the  earth’s  gravity  field. 
Hiese  equations  can  be  conveniently  divided  into  two  groups,  control 
equations  and  trajectory  equations,  which  are  related  schematically 
as  follows: 


Input  Data 


Control 

Trajectory 

I 

Equations 

Equations 

- 

Location 
Velocity 
Specific  Force 
Attitude 
Attitude  Rate 


The  control  equations  are  the  relationships  that  specify  turn  rates 
according  to  the  user's  input  data. 

The  trajectory  equations  are  a collection  of  differential  and 

algebraic  equations  that  produce  position,  velocity,  specific  force, 

* 

attitude  and  attitude  rate  in  response  to  the  imposed  control.  They 
are,  in  short,  the  equations  of  motion  for  a body  free  to  move  in 
six  directions  in  inertial  space. 

The  trajectory  equations  are  kinematic  relationships,  i.e.  they 
deal  with  motion  in  the  abstract  without  reference  to  force  or  mass. 
Since  force/mass  concepts  are  immaterial,  PROFGEN  avoids  all  aircraft- 
specific  considerations  such  as  moment  of  inertia,  aerodynamic  force 
and  thrust  force.  It  follows  that  the  aircraft  modeled  here  is  a 
weightless  body  that  can  be  displaced  and  rotated,  without  restriction, 
to  suit  the  users,  demands . 


« 


In  the  following  development  those  equations  that  became  part 
of  the  actual  code  in  PROFGEN  have  stars  (*)  beside  their  numbers. 

» 

4.1  Coordinate  System  Descriptions  and  Relationships 

The  coordinate  systems  of  particular  interest  in  this  report  are 

the  inertial,  earth,  navigation  and  path  systems.  These  four  systems, 

or  frames,  will  be  defined  shortly  as  right-handed  orthogonal  frames. 

The  relationship  of  the  earth  and  navigation  frames  will  determine 

aircraft  location  (longitude,  latitude,  alpha)  while  that  of  the 

navigation  and  path  frames  will  determine  attitude  (roll,  pitch,  yaw). 

Location  and  attitude  data  will  be  carried  in  two  direction  cosine 

matrices  (c”  and  Cn)  that  describe  the  rotations  between  pairs  of 
e p 

coordinate  frames.  The  subsequent  portions  of  this  section  describe 
the  four  frames,  define  the  two  direction  cosine  matrices  and  delineate 
the  extraction  of  location  and  attitude  angles  from  each  of  these 
matrices. 


I'-,.,  x Frame  Descriptions 

• Inertial  frame  (i  frame:  X^,  axes) 

The  inertial  frame  has  its  origin  at  the  earth’s  center 
of  mass  and  is  non-rotating  relative  to  the  stars.  This 
frame  is  important  mainly  as  it  applies  to  the  computation 
pecific  force.  Its  relationship  to  the  earth  frame 
l is  portrayed  in  Figure  9. 

L 

& 
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• Earth  frame  (e  frame:  X , Y , Z axes) 

e*  e e 

The  earth  frame  has  its  origin  at  the  earth's  center  of 
mass  and  has  axes  fixed  in  the  earth,  Figure  9.  Axes 
Y^,  Y^,  and  Z.^  all  lie  in  the  earth's  equatorial  plane 
while  axes  Xft  and  X^  are  coincident,  passing  through  both 
poles.  The  rate  of  rotation  between  these  two  frames  is 
the  earth  sidereal  rate,  designated  ft.  WGS-72  (Reference  l) 
gives  this  value  for  ft  which  is  denoted  WEI  in  PROFGEN: 

ft  « 0 . 72921151^7  x 10~k  rad/sec 

• Navigation  frame  (n  frame:  x,  y,  z axes) 

This  locally-level  frame  has  its  origin  at  the  aircraft 
center  of  mass  with  x and  y in  a plane  tangent  to  the 
reference  ellipsoid  and  z perpendicular  to  the  ellipsoid. 

Figure  9.  (Center  of  mass  and  center  of  rotation  are 

/• 

coincident  in  this  development).  PROFGEN  solves  the 
trajectory  equations  in  the  navigation  frame.  Aircraft 

location' is  specified  relative  to  the  earth  frame  by  the 

\ 

three-tuple  (X,  <|>,  a)  where  X is  longitude,  <j>  is  geographic 
latitude  and  a is  the  navigation  frame  heading  angle, 
referred  to  variously  as  alpha,  wander  angle  or  wander 
azimuth  angle.  Figure  9 shovs  that  <{>  is  geographic  latitude, 
not  geocentric  latitude.  Thus  z is  normal  to  the  elliptical 


6 


39 


M 


surface  of  the  earth  rather  than  in  the  direction  of  the 
earth  center.  The  values  for  X,  <j>,  and  a will  be  computed 
from  the  direction  cosine  matrix  C^. 

• Path  frame  (p  frame:  x , y , z axes) 

P P P 

The  path  frame,  depicted  in  Figure  10,  has  its  origin  at 
the  aircraft  center  of  mass.  It  takes  its  name  from  the 
fact  that  the  x^-axis  follows  the  aircraft  path  by  staying 
aligned  with  the  total  velocity  vector,  V.  (V,  velocity 
with  respect  to  the  earth,  will  be  defined  precisely  in 
Section  4.2.3) 

In  general  V is  misaligned  from  the  aircraft's  longitudinal 
axis  by  an  angle  of  attack  and  a crab  angle.  In  this 
development  we  assume  these  angles  are  zero.  The  effect 
of  this  assumption  is  to  weld  the  path  frame  to  the 
aircraft's  body  thus  causing  to  pass  through  the 
aircraft  nose  and  y^  to  point  out  the  right  wing. 

points  down  in  level  flight  but  rotates  about  xp 
during  coordinated  turns  so  there  is  never  any  maneuver 
acceleration  along  yp. 

Since  path  and  body  are  coincident,  the  familiar  body  frame 

terms  of  roll,  pitch  and  yaw  will  be  borrowed  to  describe 

the  Euler  angles  between  the  path  and  navigation  frames. 

Boll,  pitch  and  yaw  are  denoted  n , n and  r\  and  are 

x y z 


S 
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Note:  Origin  of  path  frame  displaced  from  that  of  nav  frame 

only  for  clarity  of  diagram;  they  are  actually  coincident 
at  aircraft  center  of  mass . 


Figure 


Navigation  and  Path  Coordinate  Frames 
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measured  around  x , y and  z respectively.  A right  turn 
P P P 

produces  a positive  yaw  rotation,  a pitch  up  is  a positive 

pitch  rotation,  and  a clockwise  roll  (as  viewed  from  behind 

the  aircraft)  is  a positive  roll  rotation.  The  values  of 

n , n and  n will  be  computed  from  the  direction  cosine 
x y z 

matrix  Cn. 

P 


4.1.2  Frame  Relationships:  Direction  Cosines  and  iuler  Angles 

# Earth  and  Navigation  Frames 

Figure  9 presents  the  relationship  between  the  earth  and  navigation 

frames.  When  A,  <J>  and  a are  zero,  the  navigation  frame  is  directionally 

aligned  with  the  earth  frame.  Beginning  at  the  aligned  position,  the 

rotations  necessary  to  go  from  earth  to  nav  coordinates  form  the 

direction  cosine  matrix  C^.  This  matrix  is  the  ordered  product  of 

three  individual  matrices  describing  these  rotations:  an  x rotation 

of  A degrees,  a y rotation  of  4>  degrees  and  a z rotation  of  a degrees. 

Using  an  "s"  prefix  for  the  trigonometric  sine  and  a "c"  prefix  for  the 

cosine,  is 
© 


Xa 

A* 

a 

-CO ( 

o 

o 

a 

t 

C<t>  O -Af 

o t o 

A#  d *<p 


/do 
o AK  AK 
o -AX 


(7) 
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Now  if  the  elements  of  sure  identified  aB 


C£W„ 

CENn 

C£N„ 

C€Nlt 

CENu 

CENn 

C£N„ 

ZEN3l 

CENU 

then  the  individual  elements  are 


CE N„  * AM* 
cenv  * “amol  am4> 

CENtl  - Am  4 

CEN,i  a Am*  AMX  /■  , codOC  Am4>  am  X, 
CENxt  * AM*  *64  X - am  oc  Am  <p  Am  X 

cENn  * * am  A 

c£W/j  “ of  <4«c  A.  - ^ar  ^ A 

C£A/*3  = /C^of  AM  A + AMOt  AM  $ *09%. 
CEN&  ' *04  4>  <04  X 


* Coded  for  implementation  in  PROFGEN. 


1*3 


5SB9HP 


To  extract  latitude,  longitude  and  alpha  from  the  elements  of  C^, 
the  following  calculations  are  made 


4>  - ( cENti)  , <f>  e[-%'  * %]  vo) 

A * &’(-  <M)t /cf^u)  , A + <"> * 

Oi m £*'(- C£N„ ) (,t* 


where  the  initial  values  for  <£,  X and  a are 

<p  a LATO 
X = LONO 
a = ALFAO 

The  FORTRAN  functions  SIN  ( • ) and  ATAN2  (*,*•)  were  used  to  implement 

(10) ,  (ll)  and  (12)  because  their  range  agrees  with  that  desired  for 
<f>,  X and  a.  An  important  aspect  of  the  computation  for  X in  Equation 

(11)  is  that  £ [-tt/2,  tt/2],  which  means  cos  (tj>)  is  always  positive, 
which  in  turn  makes  the  sign  of  CEN^  and  CEN^  depend  solely  on  X, 
which  removes  any  doubt  as  to  the  quadrant  where  X lies.  A similar 
statement  applies  to  a as  computed  in  (12). 


individual  elements  are 

i 

* cosri„  * cosn 

Z 'y 

',1 

■3 

= -sinn  * cost) 

z 'y 

] 

°ps31 

= sinn 

y 

1 

= cosn„  • Sinn  * 

z 'y 

sinn„  - sinn 

X ‘z 

• cosn 

X 

M 

1 

■« 

OT22 

® -sinn.  * sinn  * 
z y 

sinn  - cosn 
x 'z 

* cosn 

X 

OS) 

®»32 

* -cosn  • sinn 

y x 

-f 

i 

j 

^13 

- cosn„  * sinn  • 

z 'y 

co sn„  + sinn 
x z 

* sinn 

X 

’ 

I 

; 

CM23 

= co8Rz  * sinnx  - 

ainn„  * sinTV 

Z «/ 

• cosnx 

4 

j 

CW33 

* -cosn r • cosn 
y x 

a 

% 
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Roll,  pitch  and  yaw  are  extracted  from  the  elements  of  Cn  as  follows 

P 


%**■&*(-  m3t  /-cp»„)  , ^ * [-rr,trr] 

» JU*  (CPH„)  ( *[-%/%]  (,»* 

% *'&M  '(vk/cpn.)  , ^ ,fc  tw]  03)* 


where  the  initial  values  are 


n = ppitcho 

y 

n = ALFAO  + PHEADO 

Z 

Again  SIN  ( • ) and  ATAN2  ( • , • ) were  used  to  implement  (l6),  (17) 
and  (l8).  As  with  A and  a,  the  key  to  the  computations  in  (16)  and 
(l8)  lies  in  the  fact  that  n has  a restricted  range  which  makes  its 

y 

cosine  always  positive.  The  relationship  between  a,  n and  ip  (heading) 

z 

is  illustrated  in  Figure  11. 


4.2  Trajectory  Equations 

Sections  4.2.1,  4*2.2  and  4.2.3  will  develop  first  order 
differential  equations  to  describe  the  motion  of  a body  in  six 
degrees  of  freedom.  Section  4.2.4  defines  the  states  of  the 
state  vector,  x.  The  companion  algebraic  relationships  for 
specific  force,  attitude  rates  and  plumb-bob  gravity  will  be 
developed  in  Section  4.2.5. 


4.2.1  Direction  Cosine  Rates:  Location  and  Attitude 

At  least  three  methods  are  available  for  keeping  track  of  the 
rotation  angles  between  frames,  including  direct  integration  of  the 
Euler  angle  rates,  propagation  of  four  quaterion  parameters  representing 
a complete  direction  cosine  matrix  (Reference  5),  and  propagation  of 
the  direction  cosine  matrix.  The  last  approach  was  chosen  for 
PROFGEN  because  of  its  simplicity  and  versatility.  This  section 
derives  a general  expression  for  the  direction  cosine  rate  and 
then  displays  the  result  in  notation  appropriate  to  and  C^. 

For  any  two  frames,  a and  b,  the  Theorm  of  Coriolis  can  be 
written  for  any  vector  u as 


a ; 
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This  equation  is  in  "physical  vector"  form.  It  states  that  the 
time  rate  of  change  of  u,  as  observed  in  the  a frame  (i.e.  with 
respect  to  the  a framed  equals  the  time  rate  of  change  of  u,  as 
observed  in  the  b frame,  plus  the  angular  rate  of  change  of  frame 
b with  respect  to  frame  a crossed  onto  u.  The  addition  and 
multiplication  in  (19)  are  physical-vector  addition  and  physical- 
vector  cross  multiplication.  When  (19)  is  coordinatized  in  the  a 
frame,  these  "math  vector"  relationships  follow: 


vhere  is  a "cross-matrix"  that  produces  a result  on  a math 
vector  identical  to  that  of  cross  multiplication  on  a physical 
vector.  is  defined  belov.  Continuing 


Equating  (20)  and  ( 21 ) yields 


A 


'U  W 


i 


and,  since  u is  any  vector,  it  follows  that 


Uz) 
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where 


bo. 


The  specific  notation  chosen  to  implement  (22)  and  (24)  for 


Cn  and  Cn  is  shown  helow: 
e p 


0 "ft  fV 

(?3  ° "P* 

-ft  e*  ° . 


(**> 


A *> 
- 


C2*) 
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4> 


o 

*9 

O 

-*>* 

-*>4 

&■ >x 

O 

= *>, 
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For  writing  convenience.  pn  and  wn  will  be  referred  to  hereafter 
as  £ and  w.  In  (25)  and  (27)  we  have  expressions  for  keeping  track 
of  location  and  attitude  provided  £ and  w can  be  computed.  Sections 
1+.2.2  and  ^.2.3  deal  with  £.  The  computation  for  w will  be  given 
in  Section  k.3  where  turning  rates  are  discussed. 


U.2.2  Angular  Rate  - Nav  Frame  v.r.t.  Earth  Frame 

£ is  the  angular  rate  of  the  nav  frame  with  respect  to  the  earth 
frame.  The  fact  that  the  x and  y axes  of  the  nav  frame  remain  tangent 
to  the  earth  will  be  used  to  derive  expressions  for  and  py.  pg  will 
be  determined  by  the  users  choice  of  azimuth-angle  mechanization. 

Consider  the  geometry  of  Figure  12,  a section  of  the  earth 
ellipsoid,  where  and  V^.  denote  North  and  West  velocity  components. 
The  North-West-Up  (N-W-U)  frame  differs  from  the  nav  frame  only  by 
the  rotation  a.  If  V is  earth  frame  velocity,  its  navigation  frame 
components  are  denoted 


,»»  A 


(*9) 


Then  from  Figure  11,  V^,  V^,  are  given  by 


/Crtoc  o 


- o 
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The  angular  rates  required  to  keep  the  N*W«U  frame  level  over  the 
earth  ellipsoid  are  deduced  from  Figure  12  as 


~ Vw 

Rp  + h 


where  h is  altitude  above  the  ellipsoid,  R^  is  the  radius  of  curvature 

of  an  ellipsoid  meridian  line  and  R is  the  radius  of  curvature  of  the 

P 

ellipsoid  in  a plane  through  the  normal  and  at  right  angles  to  the 

meridian,  (it  can  be  shown  that  R is  the  distance  oc  where  c lies 

P 

on  the  polar  axis.)  R and  R vary  with  <P  according  to  the  following 

m p 

equations  (Ref,  2,  pp  168-170): 


R 


( | - eWo)1'* 


Es 


* 

(34) 


where 


R2_b2 

eccentricity2  - e"  ■ 0.00669** 31 7778  (WGS-72  data) 

*.2 

semimajor  earth  axis  » 209256**0  feet  (WGS-72) 
semiminar  earth  Mis  m 20855**8l  feet  (WGS-72) 


PN  and  lie  in  the  x-y  plane  of  the  nav  frame  and  can  be  resolved 
into  components  along  x and  y as  follows: 


fit  * ft. 

* - fa  A**  OC  + fa  A 

ft  * Pup 


(3S) 


The  general  relationship  between  p and  a can  be  deduced  from  the 

z 

geometry  of  Figure  12  as 


A /•  oc 


(*) 


The  value  for  a depends  on  the  azimuth  angle  mechanization  (LLMECH) 

desired  by  the  user.  The  various  choices  and  the  resulting  p values 

z 

are  tabulated  in  Table  Z • 
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LLMECH 

Name 

* 

a 

p ‘ ' 

z 

» 

1 

Alpha  Wander 

-A  sin  <p 

0 

2 

Constant  Alpha 

0 

• 

A sin  <p 

3 

Unipolar 

■ji  * 

• 

A (sin  <J>  - J) 

k 

Free  Azimuth 

-(«  + A)  sin  <f>n 

-ft  sin  <j> 

t J = sign  (<|>)  if  n * 0.729211 51^7  x 10^  rad/sec 


Table  Z - Azimuth  Angle  Mechanization  Schemes 


Figure  13,  a section  of  the  earth,  is  drawn  so  is  perpendicular 

to  the  paper  at  the  indicated  point.  (Note  again  that  R terminates 

P 

on  the  polar  axis.)  Examination  of  this  figure  shows  that  the 
equation  for  A is 


A 


-Vw 

(Rf»  + h)  /w  ^ 


p , and  p . as  derived  in  this  section,  depend  on  a,  4,  V , V and 
x y z x y 

Vz.  a and  <$>  can  be  obtained  from  c”  using  Equations  (10)  and  (12)  while 

expressions  for  V , V and  V„  will  be  derived  in  the  next  section, 
x y z 
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Figure  13  - Geometry  for  Deriving  X 


Equatorial  Axis 


1*.2.3  Velocity  v.r.t.  Earth 

Referring  now  to  Figure  9,  the  vector  R connects  the  earth's 
center  with  the  airciaft  location  at  all  times.  By  definition  of  V 
and  by  Coriolis'  Law 


(36) 


<4  )L  / _ H / J—  »v  (J9) 


Coordinatize  (39)  in  the  nav  frame  and  use  (22)  and  (28)  to  produce 
the  following  equivalent  expressions  in  math-vector  form: 
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Substituting  (h2)  in  (40)  and  writing  out  the  individual  equations 
yields 

V*  s CPN„  Vt  " 60 ^ V^  + ^3^ 

V*j  a C P N^i  Vt  + ^ Vx  * u> * N/j  f fj)* 

V^  = CPNji  Vt  “ Vx  1* 


V^  will  t>e  recognized  as  PACC,  the  path  acceleration  needed  to  alter 
the  magnitude  of  V. 

VT  = PACC  ft  /~c*  (44)* 

In  (^3)  we  have  a differential  equation  for  earth  frame  velocity  that 

T 

depends  only  on  factors  already  specified  save  for  w = (wx  0)^.  u>^)  . 

To  repeat,  w will  be  derived  in  Section  b.3 

U . 2 . h State  Vector 

PROFGEN  carries  a state  vector,  x,  'containing  23  states  in  a 23 


»*▼***«$''’*»  i.rn'i’1*' 


The  appropriate  differential  equations  for  the  elements  of  x are 

Equation  (U3)  for  the  velocity  components  V , V , V ; Equation  (M) 

x y z 

for  the  total  velocity  V_;  Equation  (25)  for  attitude  data  in  Cn; 

X P 

Equation  (27)  for  the  location  data  in  and  this  differential 
equation  for  altitude,  h; 

b = \£  (45')* 


k.2. 5 Other  Trajectory  Relationships 

The  following  three  topics  are  discussed  now  to  conclude  the 
derivation  of  the  trajectory  equations: 

a.  Specific  Force 

b . Attitude  Rates 

c . Gravity  Model 

Topic  c supports  topic  a.  Topics  a and  b are  important  only  insofar 
as  they  provide  a way  to  compute  specific  force  and  attitude  rate  for 
PROFGEN  output.  Specific  force  and  attitude  rate  are  algebraic 
expressions  not  required  during  state  vector  propagation;  therefore, 
in  some  sense,  these  equations  lie  outside  the  mainstream  of 
PROFGEN' s calculations . 


y. 


ft 

v. 
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a.  Specific  Force 

Specific  force,  F,  is  the  acceleration  that  a velocity  meter 
(accelerometer)  aboard  the  aircraft  would  detect.  Specific  force 
is  the  total  inertial  acceleration  minus  the  mass-attraction 
gravitational  acceleration;  i.e.  specific  force  is  the  second  rate 
of  change  of  R_  as  viewed  by  an  observer  fixed  in  inertial  space, 
minus  mass-attraction  gravity,  The  physical  vector  equation 

for  this  (see  Reference  3,  p.  121),  where  + and  - are  physical 
vector  operations,  is 


X 


The  navigation  frame  rotates  at  rate  T (T  = £ + fj)  with  respect  to 
the  inertial  frame.  Differentiating  (1*7) » substituting  the  result  in 
(46) , and  continuing  with  the  expansion  gives 


m 41 


'l  + -* 


= 41/  + r*y  + d.  (Q  *?)/  - §„ 

= dlj  + y * JjT  (4 X4)J  +-Q-*(Q*£)  -<Sm 

= jj/  * (p  ■'■-&)*¥  + jQ*V 


'afl+ (*+*)*  v-i 


(+8) 


63 


where  we  have  used  the  fact  dfl/dt|e=0  and  where 


? 4 Sm  “ -£•  * (-O-'B)  <W  1 

■ 1 

i 

The  vector  £ is  the  usual  plumb-bob  gravity  composed  of  both  mass  | 

attraction  and  earth  rotation  components.  The  vector  £ points  \ 

, \ 

downward.  Recalling  from  Section  4.2.3  that  | 

■ ' .1 

■ i 


we  can  componentize  (48)  in  the  nav  frame  (subscripts  x,  y,  z)  as 
follows 


'TFI*'-?  ' '■ ; V 7vT*f  ,^-T< ; ,:w^ 


where 


= C 
'-e 


cf*, 

C£N9k 


C£Nit 


(SO 


_fl  = 0.729Z//5/47  • /O  J*A/ama 


Gravity  (g^,  g^,  gn)  will  be  discussed  presently.  and  ft 


x'  y 


are  the  projections  onto  the  nav  frame  of  the  angular  velocity 
between  the  earth  and  inertial  frames ; these  quantities  are  not 
related  to  0)^,  o)y  and  ojz  except  that  both  share  the  greek  letter 
"omega",  upper  and  lower  case.  All  other  quantities  needed  for 
computing  (50)  have  already  been  discussed, 
b.  Attitude  Rates 

This  section  derives  an  expression  for  each  Euler  angle  rate 
• • • 

(o  » 0„)  as  a function  of  the  commanded  turning  rates,  w u)  and 

x y z x y 

03  . Recall  these  formulas  from  (15)  and  (27): 
z 


CPK 


CPN,. 


CPN„ 


C«?«) 


(32  i) 


(S2c) 
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CPS)a  m ~ /COO  ft *j  <03  ft). 


CPS, 


(S2tt) 


CPN„  -h  to-  cPNj,  6v«) 


CPS 3!  =- 


CPS,,  + to*  CPNt 


(S2i) 


CPS, 


-to-  CPN,3  + ti)y  cpsta 


(&a) 


Differentiate  (52c)  to  get 


cpn  = (cos  n ) n 
31  y y 


and  equate  this  to  (52f)  to  get 


-Wy  cpNii  + u)x  cpn21  = (cos  ny)  ny 


-u  (cos  n cos  n ) + w (-sin  n cos  n ) = (cos  n ) n 

Jr  Z,  Jr  X Z y V 


y y 


Assume  now  that  pitch  is  not  ±90°  and  cancel  cos  n to  yield 

y 


(S3) 


(<aofl-  P‘6) 


, • f*  ?■- 
. $ 

• >*•  'Hi 


. .rr j^ffc toifc i +<■  aattarzLUd*ic* — 


which  is  the  desired  expression  for  n . Similar  manipulations  of 

y 

(52)  produced  the  following  expressions  for  n and  n : 

X z 

C$4) 

(tmUjj*-  <0 

C$S) 

(coo  fjy 

PROFGEN  does  not  attempt  to  make  attitude  rate  calculations  when 
cos  n =0.  It  simply  prints  a warning  message  and  goes  on  (see 

y 

Section  3.3). 

c.  Gravity  Model 

Throughout  this  report  the  ellipticity  of  the  earth  has  been  accounted 

for  while  higher  order  effects  and  local  geoid  perturbations  have  been 

neglected.  The  purpose  here  is  to  derive  equations  for  gx»  g^  and 

that  are  consistent  with  this  philosophy  for  modeling  the  earth.  The 

normal  component  g will  be  tackled  first  following  the  approach 
z 

beginning  on  page  78  of  Reference  4. 

■ Derivation  for  Normal  Gravity,  g 

z 

Define  y as  gravity  normal  to  the  ellipsoid  at  altitude  zero. 

Then  for  an  altitude  h above  the  ellipsoid,  g at  this  altitude  can  be 

z 

expanded  in  a MacLaurin  series  of  terms  in  h: 


$ 

i 


1 


, ;j 


6? 


b - b(+>h) 


% 


A 


<r 


(SS.'i 


The  first  partial  is  given  by  Bran’s  formula  (Reference  1+,  Equation 
2-79)  which  is  based  on  an  ellipsoidal  earth  model: 


(S7) 


where  R and  R are  the  principle  radii  of  curvature  defined  by  (33) 
m p 


~t0-  T 


and  ( 3^ ) • 

Taking  reciprocals  and 

>/£ 

/ 

( <5^444  $ ) 

/ 

_ (/- 

» • I \ 

& — * •*) 


Truncating  these  equations,  adding  them,  and  dropping  higher  order  terms, 
produces  the  following  result 


_/ 


-L  (z  + e'-z  <£***%) 


CSV) 
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HI)*!  (»"•: 


If  y is  the  value  of  y at  the  equator  at  h*o,  the  first  order 
© 

2 

relationship  between  y and  fj  is  Jl1  = my  /R  where  m is  0.003^^9783. 

© © © 

Substituting  this  and  (58)  in  (57)*  and  simplifying,  yields 


J^j-  * ~ r V) 

2 2 

The  second  ,■  . tive  3 y3h  may  be  taken  from  the  spherical 

approximation  obtained  when  earth  flattening  is  neglected  entirely. 
Then  according  to  Newton's  law  of  mass  attraction 


(si) 


wh 


where  M is  earth's  mass  and  k is  the  uni  re r sal  gravitational  constant. 

3X  _ il  _ _ 2kM 

3h  3R  Rs 

e e 

2 2 

3 Y _ 3 Y „ 6kM 
2 2 H 

3h  3R  E 
e e 

so  that 

3h  R 

e 


CM 


6S 


Combining  (59)  and  (60)  with  (56)  produces  the  desired  approximate 
equation  for  normal  gravity. 


where  y is  gravity  au  the  ellipsoid  surface  which  is  given  in 
Reference  1,  page  22,  as 

t * de  (/+  ° oost78114  & 4-  e.oooot34i/*ui*<fi)  (U) 

<fe  * - 32. 08  77057  ft/**  (*») 

Combining  (62)  and  (63)  with  (6l),and  evaluating  all  constants, 

produced  this  final  expression  for  g : 

z 


Cjj  = - [si.0877oSr+*/69J9cr/M*p  +aooo7St8/o^%J* 

*[t°-fau27e-t-&<ortertjm4,)/,  +crsae-/s  A‘J 


■ Derivation  for  Level  Gravity,  gx  and  gy 

At  first  it  is  somewhat  surprising  to  realize  that  plumb-bob 
gravity  has  a level  component.  Such  component  arises  because  level 
surfaces  at  different  altitudes  (but  same  latitude)  are  not  parallel. 
This  fact  is  evident  when  one  considers  these  two  extremes:  at 


EH 
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h«o  the  level  surface  is  the  ellipsoid  and  gravity  points  along 
<J>;  at  the  same  latitude  but  elevated  to  h*«  , gravity  points  at 
earths  center  of  mass.  Between  these  extremes  the  difference  in 
slope  of  the  two  gravity  vectors  is  the  difference  between 
geographic  and  geocentric  latitude. 

Another  way  to  view  the  level  gravity  phenomenon  is  through 
the  curvature  of  the  normal  plumb  line  as  illustrated  in  Figure  14. 
Curvature  is  zero  in  the  east-west  direction  owing  to  the  rotational 
symmetry  of  the  ellipsoid  of  revolution.  Thus  level  gravity  is 
entirely  a north-south  acceleration. 

From  Figure  l4,  observe  the  following  relationship 


dh  * A elfi 


(*S) 


The  plumb  line's  radius  of  curvature,  r,  is  given  by  (2-22a)  in 


Reference  4: 


-L.fll 
h M 


fct) 


where  d is  distance  along  a north-south  direction.  Combining  (66) 
with  (65)  and  rearranging 


( ft  ^ 


(61) 
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The  change  in  plumb  line  direction,  8,  between  h=o  and  h=h  is 


f -+-■£&  M 

Jo  <?d 


fc*) 


An  approximate  relationship  for  d is  d = Rg<J>.  Then 


i|»=.  i$.  £* 

Id  J<t>  -?d 


Thus  (63)  becomes' 


r 


&.-L 

<?p  e* 


ft  *+ 


(69) 


To  obtain  a closed  form  expression  for  (69),  simplify  g as  follows: 

z 


$f  ~ -fi-te+e+m  -2&-*"%)h  +]jji  h ] (a,) 

£ t[/-zh/ej 

* (/jf  44*  <t>  / ft  4**  %)  ft  ~ *A/*k]  using  (62) 


Then 


' ill 
//  ^ 


~Z~  2f,4m4><eo4fy  = 


<4 


Substitute  this  in  '69)  and  integrate 


Amt jc*o4>  dh 


ZJ;  A*n  4>sC#t4> 


h 


(?o) 


3 is  the  tilt  angle  through  which  the  gravity  vector  tips  over  as 
altitude  increases.  Projecting  the  magnitude  of  gravity  (approxi- 
mated here  as  ly  |)  through  8 and  onto  the  level  surface  gives  for 
e 

gn  (g  north) 


- /63  * /O  (j)  A*t4><cc>a<t>) 


( 7 0 


Now  rotate  g^  through  a to  obtain  g x and  g^ 


n 


* & . 

n J*t<p£*)4hGiX0(  (72) 

a ^ /)  AttiP  £&4>MX0(  (/J1) 


7b 


. rnzaSl ■i.v.i.-i.  . Jj  II  ■!» 


which  may  be  stated  in  terms  of  the  elements  of  C as 

•»  , 

a = -/.  63-/0  h C£N},  C£N„ 

& = -1.63  ■'/o'* h C£N3,  C£N„ 


(74) 


(75)' 


This  derivation  of  level  gravity  was  based  on  material  in 

Section  5-6  of  Reference  4 where  it  is  pointed  out  that  the  effect 

\ 

of  topographic  irregularities  on  the  curvature  of  the  plumb  line 
often  overwhelms  the  value  from  equation  (7l) • In  high  mountains 
the  actual  deflection  could  be  10  times  greater  so  the  limitations 
of  (71)  are  apparent. 

It  • 3 Path  to  Nav  Rotation  Rates  and  Control  Equations 

The  relationships  derived  here  for  u will  produce  turning  rates 

commensurate  with  the  input  data  and  with  the  restriction  that  level- 

plane  turns  be  coordinated.  In  addition,  equations  for  controlling 

the  application  of  w will  be  derived.  This  control  will  usually  take 

the  form  of  a switch  to  turn  w on  or  off  at  a critical  event  time. 

The  control  equation  will  compute  the  event  time;  e.g.  the  time  at 

which  t}  should  be  disabled  in  a vertical  turn  to  make  An  = PITCH 

y y 

This  section  evolved  from  the  work  in  Section  3 of  Reference  6. 
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% 

rvS 

Li 


— ■ m 


HTsrTta?>? 


As  a preface,  we  list  some  basic  kinematic  equations  for  the 

illustration  below  where  S is  arc  length,  V is  speed  tangent  to  the 

path,  r is  radius  of  curvature  and  a is  acceleration  normal  to  the 

n 

curved  path: 


flight  path 


^1 

~T" de 

* 


Vs 


dS  * A * d0 


C76) 


Combining  these  equations  produces  this  relation  for  angular  rate 


<Je  _ an 

dt  ” v 


h.3.1  A General  Expression  for  w 

Now  recall  equations  (13)  and  (28)  defining  and  ok 

*['  • °i  i4f>}  °i  h • 4i  r* 


= T„c 


|o  -/  o cffj  © O too  tiff  (W 

oo-l  © o I * tlf*  Jo 


o o 
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A * A / 

£>**&/>»  « / ^ 


(ze) 


where  T. Q„ , T , T and  T are  introduced  here  for  reasons  that  will  be 
180  * z * y x 


apparent  shortly. 


Each  Euler  angle,  rix»  Hy.  and  nz»  has  an  associated  rate,  nx»  Hy 

and  nz.  For  a given  pjath  to  nay  orientation,  the  vector  associated 

with  n is  directed  along  x . If  the  given  path  frame  is  rotated  so 

that  roll  is  zero  (n  =o),  the  vector  associated  with  r\  is  directed 

x y 

along  the  new  (nv rolled)  y axis.  When  the  new  path  frame  is  rotated 

P 

* 

again  to  remove  uitch  (n  =o) , the  vector  associated  with  n is 

y z 

directed  along  the  new  (unrolled  and  unpitched)  axis.  Note  that 
these  three  vectors  are  not  mutually  orthogonal.  When  transformed 
into  the  nav  frame  and  added  vectorially,  they  give  the  entire  path 
to  nav  rotation  velocity.  Thus 


-WMTW  nn~\ 


j 


This  may  be  simplified  using  r\  - a + \p  (Figure  ll)  and  the  definitions 

z 

of  W TZ  and  CP  ln  (78): 


* cf 


Equation  (80)  is  the  most  general  relation  for  angular  rate  between 
the  path  and  nav  frames.  It  will  simplify  considerably  depending 
on  (l)  the  type  of  maneuver  (2)  the  nominal  path  (great  circle  or 
rhumb  line)  over  which  that  maneuver  is  superimposed,  and  (3)  ot  which 
is  given  in  Table  2 as  a function  of  the  nav  frame  mechanization 
choice.  In  the  following  four  subsections  it  is  assumed  that  the 
reader  is  familiar  with  Section  3.^. 


lt.3.2  Vertical  Turn 


a.  m Equation 

Since  the  aircraft's  wings  remain  level  in  a vertical  turn,  T ■ I 
and  = 0.  Thus  (80)  becomes 


-oc -iy 


78 


= ft,  * 


-OL  - <J> 


(Sf) 


where  n is  given  by  (77)  as 

y 


L - **- 

* Vrft) 

VTft)  ■ VT  (ti)  + (t-ti)  Vt 


(«3) 


a„  » race  • (pitch) 


cat) 


Note  that  TACC  is  positive  and  in  units  of  ft/sec  . A vertical  turn 
will  have  a slight  heading  rate  if  the  aircraft  is  following  a great 


circle  path  so 


• • A 

V * Vn  m 


rhumb  line 


CBS) 


t great  circie  (Section  4.3.6) 


b.  Control  Derivation 


Equation  (82)  can  be  integrated  to  yield  change  in  ri  over  the  interval 

y 

(ti>  t).  In  the  case  where  VT  varies  linearly  with  time  ( VT  = PACC  ^ 0), 
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I'i.VpI'HS fOKtt  f 


t 


(r)*tr 


t 

An 

Vr  (t{)  +(T- 


In  the  case  where  V^,  is  constant 


A%(i) 


4*-(t-tc)  ,VT 

fc  1 w 


Equations  (86)  and  (87)  may  be  inverted  to  compute  a time,  t » 
when  An  (TDONE)  = j PITCH  I : 


80 


SunlW.^^e. 


/ , Vr*o 

' (Bti 
- O (SI) 

TDONE, 

> VT  # o 

(e*f 

I ft  mo 


4.3.3  Horizontal  Turn 


a.  w Equation 

Since  the  aircraft  does  not  pitch  in  a horizontal  turn,  n is  zero 

y 

and  (80)  becomes 


* C,\°\  + 


- V 


CS9) 


where  n behaves  as  pictured  in  Figure  5.  (ii  is  either  on  or  off. 

• . • • 

When  on,  nx  = * ROLRATE.)  ip  is  the  sura  of  the  nominal  path 

♦ 

contribution  from  (85),  and  the  maneuver  contribution  due  to 


TACC: 


P = + P H 


, rhumb  line 


V*  * V*.  , great  circle 


# Coordinated  Turn  Requirement 


During  the  turn  the  normal  acceleration,  an(t),  progresses  from 

p 

zero  to  a peak  - a flat  peak  has  a magnitude  of  TACC  ft/sec  - and 

back  to  zero  (see  Figure  6).  This  progression  occurs  because  a (t) 

n 
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must  "follow"  n (t)  to  satisfy  the  1 quirement  for  coordinated  turns. 
This  requirement  manifests  itself  in  this  way: 

a„(t)  = 3z  z &*[&(*)]  (w) 

9 

Equation  (92)  shows  the  aircraft  will  turn  only  if  its  wings 

are  banked.  The  genesis  for  (92)  is  provided  in  Figure  15,  a nose-on 

view  of  the  aircraft  in  a right  turn  with  pitch  zero  (n  =0). 

y 

horizon 


Figure  15  - Balancing  Accelerations  in  a Coordinated  Turn 


The  vector  sum  of  32.2  and  an  must  act  perpendicular  to  the  wings  in 
order  to  implement  the  coordinated  turn.  Thus 

a„  -32  2 (93) 

When  pitch  is  nonzero.  Figure  15  is  altered  by  making  the  downward 
component  of  gravity  32.2* cos  n instead  of  32.2.  Equation  (92)  then 
follows  immediately. 
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# ipy  Equation 

Defining  V^  as  the  level-plane  component  of  total  speed 
(VL  * cos  riy),  we  may  plug  (92)  in  (77)  to  get  the  maneuver 
turning  rate: 

• _ 3l.l 

” ‘ vt  ( t ) 


31X 

VtM 

_ 3t  .t  ■&»  [#,  m] 

VT(t<  ) + (t-t.)  VT 


(If) 


(If) 


b.  Control  Derivation 


Examination  of  (89)  and  (9*0  shows  that  roll  and  roll  rate  must  be 
known  before  w can  be  computed.  Their  determination  rests  on 
choosing  the  appropriate  roll  history  from  Figure  5 and  then  on 
computing  TOFF*  TON  and  TDONE.  The  logic  and  calculations  for 
accomplishing  this  are  contained  in  PROFGEN  subroutines  TSETUP2 
and  YAWCHG  and  are  outlined  below. 
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• Computation 


The  most  general  roll  history  is  pictured  in  Figure  16.  (This 
roll  history  is  for  a right  turn.  A left  turn  would  be  the  negative 
of  Figure  16). 


Figure  16  - Roll  Angle  History  (Case  A) 


We  wish  to  compute  the  change  in  heading,  A&„  that  would  occur  if 

Figure  16  was  the  roll  history.  For  this  purpose  we  may  assume  t^, 

t2  and  t^  are  time  increments  measured  from  a ^ of  zero.  Set  up  the 
integral  of  (9h)  as  follows: 


Mm  s 3i. 


81* 


Prom  (92) 


JL  $ 


j &H/VA++  I 


- TACC 


%#*  31.Z  3ZZ<c*ifL 


(V) 


Recalling  (83)  for  VT.(t),  the  middle  integral  in  (96)  is 


u 

MrCr) 


3lt£*4j  W / 


frt)w7  y 

Vrft.)/ 


Vr*o 


(??8) 


***  ,Vt-o 

3t.2  m»iL  St, 


The  first  and  third  integrals  in  (96)  are  not  closed-form  integrable  unless  VT*0. 
Satisfactory  approximations  have  been  obtained  for  them  by  substituting 
VT,  average  speed,  for  V^(t)  as  follows: 


f 'it 

Jo  \lr 


&*tU(r) 

Vt(T) 


~ f X-  4 (r)Jr 

vn  Jo 

= -=f~  -&«(nx  r)Jr 

Vri  /o 

~.lo  [■*»  (lit  t)J 


Vr/ 


where  = ROLRATE  and  VT1,  average  speed  in  (t^  t ),  is 

vr(  - Vrfej  _ ^ (ie)  + tVr  &oe 


Similarly 


f^A. &&)  , ~ 

Jt,  MT (r)  * 


(to/) 


\ Itj  = ~'U)  2 VT(ij)  - Vr  tti)  +(U  + t)^  0o,) 


Inserting  (98)  - (101)  in  (96)  and  simplifying  gives 


A<//„  [■<*> (ri.t,)]  + 

**  !**>+ 

I r %T*m*t/Vr  f 

J , TACC  J«U+  Vr(t)  J S.  . 

\ + — — r •>  VT*o 


Vr 


“ + tacc  Cti~ t,)  Vr 

Vr  ^ 


Since  n.irTnPY  and  nx  ( =ROLRATE ) are  known,  t is 


An  ( TACc/3i.Zuoffu) 


• Reasoning  on  ^max 


PROFGEN  determines  if  the  maneuver  can  be  completed  (HEAD  reached) 

by  seeing  how  far  the  aircraft  would  turn  if  turn  acceleration  was  left 

on  for  the  entire  segment,  t^  to  t . Equation  (103)  is  used  for  this 

purpose  where  t^  is  obtained  from  (10*0  and  tg  is  placed  t^  seconds 

short  of  t,  = t„.  (if  2 t,  exceeds  SEGLNT,  t,  is  set  to  SEGLNT  t 2.) 
d f 1 l 

PROFGEN  solves  for  Ail)..  in  subroutine  YAWCHG  using  (103). 

Mmax 

If  exceeds  |HEAD|  , the  turn  can  be  completed  and  either 

Case  A or  3 of  Figure  5 is  appropriate.  Having  decided  A or  B (not 
C or  D) , the  problem  becomes  determination  of  t^  and  tg 
(In  Case  B,  t^  = tg.)  The  following  paragraphs  will  der^.  .e  equations 
for  t^  and  tg  for  both  Case  A and  Case  B. 

If  AiJ^^  falls  short  of  |HEADj  , the  turn  cannot  be  completed  and 
either  Case  C or  D of  Figure  5 is  appropriate.  For  Cases  C and  D, 
the  determination  of  t^  and  tg  is  trivial. 


• Case  A Roll  History 

The  roll  history  shown  in  Figure  16  is  identifiable  as  Case  A 
from  Figure  5.  Setting  A = |HEAD|  and  obtaining  t^  from  (10U), 
everything  is  known  in  (103)  except  tg,  the  time  at  which  roll-out 


should  begin.  Unfortunately,  (103)  cannot  be  easily  inverted  for  t^* 
This  difficulty  was  overcome  by  ridding  (103)  of  its  "in"  function 
which  was  accomplished  by  approximating  the  middle  integral  of  (96) 
just  as  the  first  and  last  integrals  in  (?6)  were  approximated  earlier. 
Thus  (98)  becomes 


(tn~  tj 

322  xvorfp  Vrz 


t/oS) 


where 

Vt2  = V(t,)  + 


Oex) 


Now  replace  the  first  term  in  (103)  with  (105)  - (106),  set  = |HEAD| 
and  simplify  to  get  this  quadratic  equation  in  t^: 

t][(i)k\lr  +(i)L\Jr/b,  - TAccj  Vr 

* t<[ ' '£) 4 L Vt+zLVt- TAcebt  / TAcc t,  i/TJ 

+bt[u,  i2i  + TAcc  tj  — O Oo7? 
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in  irimr — 


32.2  too  fly  A (fl%  tt)j  ~ 

lne*o\  A&rfp 
Vtti)  + t,  <JT/Z 


Ooa  ) 


Note  that  (107)  reduces  to  a linear  equation  in  tg  if  V is  zero. 

The  coefficients  in  (10 7)  are  computed  in  TSETUP2  and  supplied 
to  QUADRT  where  tg  is  computed.  TOFF,  TON  and  TDONE  are  given  below. 

To  reference  them  to  true  time  instead  of  t.=0,  merely  add  TI  to  each  one. 


ToFF  e t,  * JLf* [tACc/ (31.2 coo tfj)] 
TO*  - ~ ^rdJLi  C/07) 

TDO/VF  xt</  ~ tz  + t, 


C/o4) 


Case  B Roll  History 


A "Case  B"  type  roll  history  is  illustrated  below. 


Roll  Into  Turn  Roll  Out  of  Turn| 


ti(Tl) 


1 1 (ton=toff) 


t . (TDONE) 
a 


- Roll  Angle 


tf (tf) 


The  following  equation  in  was  obtained  uBing  a procedure  like  that 
which  lead  to  (103): 

&4A  A k*  \h&d\[ VTt,  + yr(t)j-  o (/of) 

\ 

2 

If  In  (cos  x)  is  approximated  as  -.632x  , (109)  becomes 

t,  [-40.1  it]  + t[ilMA0\rtJ+  Vrd.  ) I HfAoj  fjyc  *0  (//O) 

The  coefficients  in  (110)  are  computed  in  TSETUP2  and  supplied  to 
QUADRT  ./here  t^  is  computed.  TOFF,  TON  and  TDONE  follow  immediately 
(see  Figure  17)  when  t^  is  known. 

k.3.^*  Sine  Maneuver 
a.  (jJ  Equation 

Since  the  aircraft  does  not  change  pitch  in  a sine  maneuver,  is  zero. 

Hence ^ in  (80)  reduces  to  a form  identical  to  that  for  a horizontal  turn: 

i 


ip  is  the  stun  of  the  nominal  path  contribution  from  (85),  and  ip^, 


the  maneuver  contribution  due  to  a (t): 

n 


= (t'M  + % 


, rhumb  line 


+ V6  , great  circle 


The  next  two  paragraphs  derive  expressions  for  ik,  and  n . Expressions 
" M x 

ford  and  ijj  are  given  in  Table  2 and  Section  4.3.6,  respectively. 


b.  \pM  Equation 


For  the  aircraft  to  fly  a sine  maneuver,  its  heading  must  vary  per 


Equation  (5): 


- f+  A A**Mrt  t i^t  < Tp/l 
- A , tyt  i Tp 


where 


A = max  heading  variation 

w = heading  oscillation  frequency 

T = 2tt/w  = period  of  one  full  oscillation 


Differentiating  (5)  twice  gives 

«V>  = i 

A Ax  (2 ms  t ) 

- A MS  Am  (iM/t) 

p 

V*  = 

I z A AiT/C#  (iMxi) 

l 

"ZAjUJ  A66  (Za tst) 

• 

c.  n Equation 

In  the  context  of  a sine  maneuver,  (77)  becomes 

_ a*(i) 

a„(t) 

^4  (t)  'C&fh  \ 

where  a^Ct)  acts  in  the  level  plane  and  is  level-plane  speed*  No* 
recall  (92),  the  coordinated  turn  requirement  relating  »n(t)  to  nx(t) 

Q-n  (t)  = 322  Axarfj  (t)]  M 


Plug  (92)  into  (113)  and  rearrange  to  get 


from  which  it  fellows  that 


32.2  Vt 

32ZX  3-  (Vrib)' 


</>M 


on* 


The  trouble  that  was  experienced  in  establishing  roll  control  in  the 
horizontal  turn  is  entirely  avoided  in  the  sine  maneuver  because  there 
is  no  need  to  compute  any  special  "event  times". 


4.3-5  Straight  Flight 


a...w  Equation 

Since  aircraft  attitude  remains  fixed  in  straight  flight,  (80) 
simplifies  to 


where 


u\ 


V 


1 

\un 

i>  = <i> n = 


O/t) 


rhumb  line 
great  circle 


(as) 


To  repeat,  a is  given  in  Table  Z . As  with  the  sine  maneuver,  no 
"events"  occur  in  straight  flight  so  (ll6)  tells  the  whole  story. 


It. 3.6  Heading  Angle  Turning  Rate  for  a Great  Circle  Path 

Figure  18  shows  the  geometry  associated  with  the  problem  of  determining 

the  rate  of  change  of  heading  along  a great  circle  route.  (The  E frame 

in  Figure  18  is  established  here  to  facilitate  this  analysis).  A great 

circle  route  lies  in  a single  plane,  Plane  I,  which  passes  through  the 

center  of  the  earth.  This  plane  is  described  by  A and  ip  where  A 

eq  eq  eq 

is  the  longitude  at  which  the  great  circle  plane,  Plane  I,  intersects 
the  equatorial  plane  and  \|»  is  the  heading  at  the  aforementioned 
intersection. 

Consider  a vehicle  at  point  P proceeding  along  a path  lying  in 

Plane  I.  The  coordinates  of  this  point  are  given  by  A-Ag^,  4>c  and  R 

where  <f>  is  the  geocentric  latitude  and  R is  the  length  of  the 
c 

geocentric  radius  vector. 

The  geocentric  heading,  ip^,  at  point  P is  given  as  the  angle 

between  the  horizontal  velocity  vector  and  a vertical  plane.  Plane  II, 

erected  at  longitude  A- A and  containing  point  P.  Thus  ^ is  the 

eq  c 

angle  between  Planes  I & II.  In  rectangular  coordinatesfXg,  , ZE} 
the  equations  for  Planes  I & II  are  respectively 


0t7) 


ye  Pjj  - o 


The  angl^  is  therefore  given  by 


r 


V>C  - tfUj  (/. -Ay)  0/9) 

The  primary  concern  here  is  with  the  geographic  heading  angle  <Ji 
rather  than  the  geocentric  heading  These  two  angles  differ  due 

to  the  deviation  angle  D between  the  local  vertical  and  the  geocentric 
position  vector,  see  Figure  18.  If  <J>  denotes  geographic  latitude ,>then 

D * <t>  - Oio) 


The  projection  of  onto  a local  level  coordinate  system  (rotation 
through  angle  D about  the  local  east  axis)  yieldB 

Ain  v>< 

1 ~ **4  %aH*D 

and 


Crt  V - 


<64  D 


On ) 


Differentiation  of  (121)  with  respect  to  tine  yields  the  desired 
quantity: 


(p*£rtP  -h 
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The  Remaining  steps  are  concerned  with  determining  U9able  expressions 
for  the  right  side  of  (123) v 

The  time  derivative  of  i|»c  is  obtained  from  equation  (119)  as 


(pt,  - A vAao  - <coo  <p* 


024) 


From  Figure  18  it  is  seen  that 


]/f  ~ fa 


Off) 


}?£  *“  fi?  AAO  sCOO  {'A  ) 024) 


Combining  equations  (UT),  (118),  (119),  (125)  and  ( 126 ) yielus 

» 2 2 • | 

M A*>4>*  + Am  <pA  ( 

which  when  substituted  into  (12^)  gives  the  simple  expression 

* A AMt  (<p-D)  < 


027) 


028) 


Also  required  is  the  inverse  solution  of  equations  (121)  and  (122) 
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AM  V <XOD 

Am1 0 + <c*0*O 


028) 


A*0  V 
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Substituting  (128),  (129)  and  ('  30)  into  (123)  yields  the  expression 
for  the  turning  rate  of  the  heading  angle 

(f)  = As**  ($‘D)  A# D 03/) 

Since  this  is  the  value  of  ip  required  to  maintain  flight  in  the  great 

circle  plane,  it  is  the  quantity  labeled  previously  as  ip  . Thus 

G 

ipG  <=>  Equation  (l3l) 

The  angle  D and  its  time  derivative  D are  given  for  the  ellip- 
soidal earth  by  the  following  relationships: 


o » 
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Equations  (131)  through  (135)  are  the  exact  relationships  for 
an  ellipsoidal  earth.  If  the  earth  had  been  assumed  spherical,  its 
eccentricity  would  have  been  zero  and  (131)  through  (135)  would 
reduce  + / 

iG  X sin  <j>  (!36) 

D = 0 

R = R + h 
e 

D = 0 

• • 

R = h 
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The  earth  model  in  PROFGEN  may  be  converted  from  an  ellipsoid  to  a 
sphere  by  simply  seating  ez=0  in  BLOCK  DATA.  However,  to  take  full 
advantage  of  the  spherical-earth  simplification  would  require  replacing 
(131)  through  (135)  with  (136)  and  revising  the  gravity  and  earth  radii 
computations . 


SECTION  V 


PROGRAM  ORGANIZATION 

Previous  sections  have  described  what  PROFGEN  does , explained 
how  to  use  it,  and  derived  equations  for  its  implementation.  This 
section  assembles  these  equations  in  a sequence  amenable  to  solution 
in  FORTRAN  code.  Flow  of  equations  and  code  are  both  presented. 

Two  principles  will  guide  us  now  (Reference  7): 

• The  most  reliable  documentation  for  any  program  is  the  code  itself. 
Therefore  our  purpose  is  not  to  describe  the  code  in  minute  detail  - 
such  a description  would  be  unreliable,  redundant,  and  probably  harder 
to  read  than  the  code  itself  - but  merely  to  show  how  large  pieces  of 
code  interact. 

• Each  subprogram  contains  comments  giving  a readable  description  of 
what  that  subprogram  is  supposed  to  do.  These  comments  form  the  core 
of  the  micro-level  documentation  and  do  not  need  to  be  repeated  here. 

Figure  19  is  a macro-level  flow  chart  emphasizing  overall 
computational  structure,  especially  with  regard  to  control  of  step 
size,  h.  The  name(s)  beside  each  block  in  Figure  19  designates  the 
subprogram( s ) where  the  action  in  that  block  occurs.  Each  of  these 
subprograms  usually  cedis  one  or  more  other  subprograms  to  complete 
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this  action  (see  Figure  21).  The  main  program  and  master  executive 

* 

is  named  PROFGEN.  The  subexecutive  for  controlling  numerical 
integration  during  each  maneuver  is  FLTPATH. 

Figure  20  is  an  expansion  of  the  integration  block  that  appears 
in  heavy  outline  in  Figure  19.  Figure  20  was  included  here  to  show 
how  the  differential  equations  in  Section  IV  actually  get  solved. 

Figure  21  is  a dependency  chart  showing  what  calls  what. 
Although  timing  relationships  are  vague,  (Figures  19  and  20  deal 
with  timing)  this  chart  is  nevertheless  useful  for  getting  a bigger 
picture  of  how  PROFGEN  fits  together.  It  was  kept  during  program 
development  to  help  assess  the  impact  of  proposed  changes. 
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Figure  19  - Macro- Level  Logic  Flov  Diagram 
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r Enter  KUTMER  > 
. from  FLTPATH  v 


Save  x and  t 
— o o 


V0=X  ( rn ) 

VlsY0t(H/*)F(T0,V0J 

Y2=V0MM/ ■■>)F(TO,YO)  ♦ < 1H/6 ) F ( TO+ H/~ , Y i ) 

Y :<■=  YC+  ( **/  ) F ( TO,  YO  1 + C F (TO  + H/  ’> , Y?) 

Y't-=YO+(u/?)fr(tT,YO)-(.?H/^)r(TntH/V,Y,Z)»-(?H)r'(TO+H/?,Y’) 

Y?  - YO*'  (H/f))  c ( f'1»  YO)  + (2H/1 ) r t y ^ ) t (M/f) ) r (t  y/,  i 

~y  (T O + H ) 


Note: 

All  computations 
except  F(tt  x) 
occur  in  KUTMER, 


Figure  20  - Numerical  Integration  of  x ■ F(t,  x)  from  tQ  to  tQ+h 
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VALDATA  1TEVUHIT  HEADEK  KMPERF  FLTPATH  PRTPOUT  BITEC  UT 
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fDISSPLA  Subprograms 
(these  sources  are  not 
listed  in  Appendix  B) 
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Figure  21  ( Continued) 


'■H 


3 


106 


APPENDIX  A 


SAMPLE  RUN  OF  PROFGEN 

The  sample  run  described,  here  was  constructed  in  seventeen 
segments  to  exercise  moat  of  PROFGEN’ s code  including  at  least  two 
segments  of  each  of  the  four  types  of  maneuvers.  Throughout  the 
sample  run  the  nominal  flight  path  is  a great  circle,  the  output 
interval  is  one  second,  and  the  integration  step-size  is  variable. 
Figures  2 and  3 are  the  PRDATA  and  PASDATA  lists  that  were  used  as 
input . 

• Printed  Output 

Figure  A-l  is  a •nortion  of  the  printed  output  from  the  sample 
run.  The  first  page  of  printed  output  consists  of  a banner 
(automatically  printed  by  the  CYBER-71*  computer)  followed  by  the 
date  and  clock  time  of  the  run.  The  second  and  third  pages  are 
listings  of  the  PRDATA  and  PASDATA  lists  as  read  from  TAPES? » the 
local  file  o.'  which  the  input  should  reside . These  listings  simply 
echo  thn  '.ata,  including  its  mistakes  if  any. 

Page  of  Figure  A-l  begins  the  printed  output  generated  during 
the  computational  portion  of  the  ran.  with  IPRNT»l.  This  output 
consists  of  a header  *nl  t list  of  variable  values  at  the  start  of 
each  segment,  followed  by  output  at  DTO  intervals  (one  second  in  this 
run)  during  the  segment.  The  list  of  variables  printed  does  not  change 
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and  the  definition  for  each  3uch  variable,  with  its  units,  is  given  in 
Table  A-l.  Pages  5,  6 and  7 of  Figure  A-l  show  output  up  through  the 
beginning  of  segment  2. 

The  last  page  of  the  sample  printc.it  contains,  in  addition  to 
output  spaced  at  DTO  intervals,  output  at  t-final  (460*5  seconds  in 
this  run)  plus  a post-run  assessment  of  the  numerical  integration 
burden.  In  this  case  5393  numerical  integration  steps  were  used  and 
F was  called  34675  times. 

• Plotted  Output 

Figures  A-2  through  A-6  are  the  plotted  output  for  the  sample 
run.  The  small  numbers  appearing  along  the  curve  in  each  figure  are 
segment  numbers  designating  approximately  where  each  new  segment 
began.  The  latitude  - longitude  plot  in  Figure  A-2  is  constructed 
with  the  latitud.  and  longitude  axes  at  the  same  s-'ale. 

• Other  Output 

TAPE3  output  was  suppressed  in  the  sample  run  by  setting  IRITE*0, 
If  TAPE3  output  had  been  specified  (unformatted  binary  records),  each 
record  would  havo  contained  the  following  list  of  variables  in  units 
of  feet,  seconds  and/or  radians:  time,  latitude,  longitude,  alpha, 

altitude,  roll,  pitch,  yaw,  velocity  components  along  nav  x,  y,  z 
and  specific  force  components  along  nav  x,  y,  z.  Subroutine  RITEOUT 
should  be  consulted  if  a more  definitive  description  of  TAPE3  output 


is  needed. 


Table  A-l 


Output  Variables 


Variable 

Unite 

Description 

..  . . . . 

TIME 

sec. 

time  (t) 

LAT 

deg. 

geographic  latitude  (<t>) 

‘ 

; 

LON 

deg. 

longitude  (X) 

K 

ALPHA 

to 

• 

■cJ 

angle  between  north  and  nav  X-axis  (oi) 

ALT 

feet 

altitude  from  ellipsoid  (h) 

ROLL 

deg. 

roll  (nx) 

PITCH 

deg. 

pitch  (n  ) , 

YAW 

deg. 

yaw  (nz) 

PSI 

deg. 

ground  heading  angle  measured  positive 

cw  from  north  (ip)  j 

DROLL 

deg/eec 

derivative  of  roll  (nx) 

; 

: 

DPITCH 

deg/ sec 

derivative  of  pitch  (n  ) 

y 

J 

DYAW 

deg/sec 

derivative  of  yaw  (n  ) I 

Z | 

VX 

ft/sec 

velocity  w.r.t.  earth  along  nav  x-axis  (V  ) 

V 

VY 

ft/sec 

velocity  w.r.t.  earth  along  nav  y-axis  (V  ) I 

v 1 

VZ 

ft/sec 

velocity  w.r.t.  earth  along  nav  z-axis  (V^)  j 

VPATH 

ft/sec 

magnitude  of  total  velocity  (V  ) 1 

A i 
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APATH 
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